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SOLID  STATE  OPTICS 


Tnis  contract  concerned  itself  with  two  major  topics 


Nonlinear  Interaction  of  Light  with  Matter 


BRIEF  SUMMARY 


rforK  on  nonlinear  interaction  of  light  with  anti- 
ferromagnetic materials,  and  thermoelastic  materials  by  tty 
ho  is  described  in  more  detail  below. 


Tne  Dirac  Theory  of  Poisson  brackets  usually  used  to  handle 
the  constraint  problem  in  electrodynamics  has  been  applied 
oy  C.  Valenti  to  handle  magnetic  materials  whose  magnetiza- 
tion vector  has  a fixed  length. 


A microscopic  approach  toward  nonlinear  optics  has  been 
started,  by  C.  Huang.  His  work  reported  here,  provides  a 
localized  oroital  estimate  for  nonlinear  optical  constants. 


Tne  paraxial  wave  optics'  method  with  intensity  and  slope  as 
the  Key  variables  avoids  the  rapid  phase  fluctuations  that 
occur  when  field  amplitudes  are  used  as  variables.  A pro- 
gram to  use  equations  for  these  variables  to  describe  an 
unstable  resonator  with  rectangular  mirrors  is  being  written 
in  collator  tion  with  Professor  W.  H.  Louisell  of  the 
University  of  Southern  California. 


Tne  spatial  distribution  of  the  temperature  rise  in  a solid 
induced  by  a Gaussian  laser  beam  of  width  w and  power  P is 
expressed  as  a one  dimensional  integral.  The  maximum  tem- 
perature rise  is  given  by  the  closed  form  expression 


max 


where  K is  the  thermal  conductivity,  N(w)  is  a 
factor  (<  1 ) that  depends  on 

W = aw  s (power  attenuation  cons tant) ( beam  width) 


March  1 0,  1 977 


r 


! 


- 3 - 


= (beam  width)/ (attenuation  length). 


The  closed  form  expression  for  the  reduction  factor  is 

J.  *2 

N(W)  = WD(^W) — e 4 Ei(^  W2) 

2 2 \ ! n 4 

wnere 

? ,°°  5 

D(x)  = exp  (-xz)  J exp  tr  dt 

0 

and 

Ei(x)  = P J'°°  dt 

-oo 

the  tabulated  functions2,  Dawson's  integral  and  the  exponen- 
tial integral  respectively.  A plot  oF  N(W)  is  given  in  Fig. 

1. 

As  an  example,  a 5p  beam  with  an  attenuation  length  i/a  - ip 
in  GaAs  has  W = 5,  N = .83  0 and  a maximum  temperature  rise  of 
(1  24)  (.83)  = 1 03°K. 

Electrodynamics  in  Thermoelastic  Die lect rics 


An  ab  initio  long  wavelength  Lagrangian  theory  of  electro- 
dynamics of  thermoelastic  dielectrics  was  developed.  It's  pur- 
pose i s to  extend  the  Lagrangian  method  developed  by  Lax  and  Nel- 
son4' » to  include  thermal  effects.  In  general,  this  is  the 
treatment  of  the  interactions  of  electromagnetic  fields,  acoustic 
waves  and  heat  conduction  with  dissipation  in  dielectrics. 


Our  approach  is  to  construct  a Lagrangian  density  which  con- 
sists of  the  field  Lagrangian,  the  interaction  Lagrangian  and  the 
matter  Lagrangian.  The  matter  Lagrangian  is  equal  to  the  Kinetic 
energy  less  the  internal  energy  U which  does  not  arise  from  clas- 
sical electromagnetic  forces.  The  internal  energy  is  a function 
of  the  finite  strain  tensor,  body  components  of  the  excitation 
fields  ana  the  entropy  density  S.  The  absolute  temperature  T is 
defined  by  T = dJ/db,  and  in  general  is  a function  of  the  posi- 
tion vectors.  The  free  energy  density  F is  then  defined  by 
F » U - TS,  so  F is  a function  of  the  finite  strain  tensor,  body 
components  of  tne  excitation  fields  and  the  absolute  temperature, 
by  the  definition  of  the  free  energy,  we  have  the  relations: 


( Aii ) 

'S 

Aa 


( ALL. 

dyj 
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wnere  are  the  components  of  finite  strain  tensor  and  are 

componuents  of  internal  excitation  fields.  Therefore,  the  Stress 
tensor  ana  the  thermal  expansion  force  acting  on  the  particles 
can  he  derived  from  the  free  energy  F.  As  far  as  the  irreversi- 
ole  processes  are  concerned,  there  is  a friction  force  which  can 
not  arise  from  the  potential  energy.  A dissipation  function  D 
can  oe  defined7 »d  by  the  rate  change  of  the  internal  entropy  pro- 
uuction.  The  dissipation  function  D is  a function  of  the  rate 
change  of  deformation,  the  rate  change  of  body  components  of  the 
excitation  fields  ana  the  temperature  gradient.  All  equations  of 
motions  for  the  fields  and  matter  with  thermal  effects  and  the 
constitutive  equation  of  Ouhamel -Neumann  Law9  can  be  derived  from 
tne  Lagrangian  and  the  dissipation  function  without  a phenomeno- 
logical formulation.  The  conservation  laws  of  momentum,  angular 
momentum  ana  energy  are  automatically  satisfied  from  the  field 
equations  ana  the  invariance  properties  of  the  free  energy  F and 
tne  dissipation  function  0.  The  equation  of  heat  conduction  cou- 
pled with  deformation  and  excitation  fields  is  also  discussed. 
In  further  work  we  will  specialize  to  linear  thermolasticity  and 
piezoelectricity  and  intend  to  apply  these  results  to  ferromag- 
netics, antiferromagnetics  and  liquid  crystals. 

Pnoton-Maunon  Interaction  in  Anisotropic  Magne tic  Crvsta Is 

Tne  coupling  between  photons  and  magnons  in  anisotropic  far- 
romagnets  ana  antiferromagnets  has  been  studied.  By  solving  the 
equations  of  magnetic  field  H from  Maxwell's  equations  and  of 
magnetization  M under  the  effects  of  constant  and  rf  fields,  we 
shall  be  able  to  obtain  the  dispersion  relation  of  the  coupled 
system  in  anisotropic  crystals.  When  the  proper  boundary  condi- 
tions are  applied,  we  can  obtain  the  coefficients  of  reflectivity 
and  transmissivity.  _Both  coefficients  depend  on  the  external 
static  magnetic  field  and  the  thickness  of  the  film.  The  pos- 
siole  application  of  using  a thin-film  anti  ferromagnetic  crystal 
as  an  optical  device  to  obtain  pure  circularly  polarized  radia- 
tion. jhy.  transmission  has  been  discussed  before  by  Zuniga  and 
bose  * . Tne  Green's  function  method  developed  oy  Lax  and  Nel- 

son4* was  applied  to  solve  the  magnetic  field  equation  in 
anisotropic  crystals.  The  results  were  used  to  obtain  the 
aosorption  or  scattering  power  which  can  be  applied  to  study  the 
spin  resonance  spectra  in  anisotropic  ferromagnets  and  antifer- 
romagneta. 

Our  treatment  starts  from  the  coupled  equations: 

H J i2) 

K2  - *2/c2  * ' 
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= r M x [H0  + H + HetC]  , (3) 

— gf  f 

where  H is  the  effective  crystal  field  and  T is  the  gyromag- 
netic  ratio.  The  abbreviated  notations  k = -iv  and  ui  = id/dt  are 
used  to  get  Eq.  (2).  The  magnetization  can  be  solved  in  terms  of 
the  magnetic  field,  M = x * ri»  using  Eq.  (3),  where  xm  is  the 
magnetic  susceptibility  tensor  which  is  a function  of  frequency, 
constant  magnetic  field  HQ  and  anisotropy  field  of  the  crystal. 
M can  be  substituted  into  Eq.  (2)  to  obtain  the  dispersion  rela- 
tion of  photon-magnon  interaction.  If  we  are  looking  for  the 
scattering  power,  E q.  (2)  can  be  rewritten  as 

o(k,u*)  * h(r,aj)  = MNL(r,i»>)  (4) 

a(k,w)  = (£)2  [ k2  I - k k ] - £m(u>)  (5) 

where  is  the  frequency  dependent  magnetic  permeability  ten- 
sor is  the  nonlinear  part  of  the  magnetization.  Eq.  (4) 

has  the  Green's  function  solution1  J 

HU,u))  = j‘  G(r  - r‘  ) * MNL(r')  dr'  e-i,ut  , (6) 


where 

a(K)  = j «*p  u*..’  «)  , 

2(k,ui)  (2«)J 


(7) 


When  the  magnetic  field  in  Eq.  (7)  has  been  solved,  it  is  easy  to 
calculate  the  electric  field  and  the  scattering  power  in  terms  of 
the  magnetic  permeability  tensors. 

Tnese  results  obtained  here  will  be  applied  to  the  ferromag- 
netic ana  antiferromagnetic  crystals  with  uniaxial  and  biaxial 
structures.  In  further  work,  the  photons,  magnons  and  photons 
interaction  in  anisotropic  magnetic  crystals  will  be  studied,  and 
tne  interaction  of  magnons  with  electric  field  will  also  be  con- 
sidered. 

Application  of  Ui rac  Constraint  Theory  to  Magnetism 

In  working  out  a consistent  macroscopic  Lagrangian  theory  of 
linear  and  nonlinear  electrodynamics  for  an  anisotropic  magnetic 
dielectric  possessing  acoustic,  ionic,  electronic  and  spin  exci- 
tations, it  is  necessary  to  identify  the  classical  coordinates 
corresponding  to  the  spin  degrees  of  freedom.  This  identifica- 
tion is  needed  in  order  to  obtain  the  Euler-Lag range  equation 
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relating  to  the  spin  motion,  and  so  that  a complete  canonical 
stress  tensor  containing  the  spin  contributions  may  be  derived. 

A classical  Lagrangian  and  Hamiltonian  theory  of  a rigid 
magnetic  continuum,  easily  generalized  to  the  elastic  case,  has 
been  formulated  which  indicates  that  the  choice  of  the  magnetiza- 
tion m(z,t)  ana  its  time  rate  of  change  m(z,t)  as  generalized 
coordinates  ana  velocities  corresponding  to  the  spin  degrees  of 
freedom,  is  a proper  one.  This  calculation  also  arouses  interest 
in  that  Dirac's  generalized  dynamics,  hitherto  used  in  connection 
with  quantum  electrodynamics  and  quantum  field  theories,  finds  an 
application  in  the  physics  of  the  continuum  solid. 


Secona  Order  Optical  Susceptibilities  of  III  - V Compounds 


The  long  wave  length  second  order  optical  susceptibilities 
of  tetrahedrally  bonded  materials  can  be  calculated  within  the 
context  of  the  pond  orbital  model  , which  is  an  empirical  tight 
Dinding  model.  By  the  symmetry  of  the  zinc  blende  structure,  the 
only  nonvanishing ^components  of  the  second  order  optical  suscep- 
tioilities  is  x^3« 


3 

The  electrons  will  gain  a potential  energy  U = £ e Ei  xi 

per  bond  when  an  external  electric  field  e is  applied  to  the 
solid,  where  x.  is  the  coordinate  measured  from  the  center  of  the 
bona  and  is  along  the  three  cubic  edges.  The  ground  state  energy 
By  of  one  bond  can  be  expanded  as 


(a) 


Bach  expanding  terms  can  be  expressed  by  the  usual  perturbation 
tneory  as 


*0)  .. 


< GlxjG  > = e Ei  XA  , 


42)  * ' ^ Ei  JSj  2- 

i,j  1 3 M 


- . <G|Xi|M>  <M | X j | G> 

Torr^J— 


(9) 

(10) 


M 


*13  > <G|X1|LXL|X,  - X,|MXM|XK|S> 

U 4k  1 “i  “*  t.K  <E<°>  - E<°>)  (E<°>  - E<°>)  ' (,1> 

ii  vj  M vj 


where  the  summations  of  L and  M are  over  all  the  exciting  states. 
Comparing  the  above  equations  with  the  total  macroscopic  elec- 
trostatic energy: 
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where  N is  the  number  of  bonds  in  unit  volume  and  § is  along  the 
body  diagonal  direction  of  the  unit  cube. 

in  terms  of  the  parameters  defined  in  the  bond-orbital 
model  , x^23  can  De  simply  expressed  as 

X123  = ^ ieix(1  ) d T ap  / E2  , (15) 


where  is  the  principal  optical  absorption  peak, 
length,  a is  the  polarity  of  the  bond  and  T is  a 
ter  for  each  row  of  the  Periodic  table. 


d is  the  bond 
fitting  pa ra me- 


in table  I,  the  values  of  x{  for  all  the  III  - V compounds 
evaluated  here  are  listed  in  column  I.  For  comparison  the  values 
obtained  from  the  dielectric  model  of  J.  C.  Phillips  and  J.  A. 
Van  Vechten  are  listed  in  column  II,  the  values  of  D.  A.  Kleinman 
using  the  same  model  are  listed  in  column  III.  In  column  IV  and 
V are  the  values  obtained  by  the  bond  charge  model  of  B.  F. 
Levine  ana  the  experimental  results.  Due  to  the  uncertainties  in 
tne  measurement  of  the  absolute  values  of  xi  “ ' experimentally  the 
relative  values  or  the  trends  may  be  more  ^important  in  these 
xinas  of  simple  calculations.  The  correct  trends  are  obtained 
here  comparing  with  the  experiment  as  well  as  the  previous 
authors.  However,  in  the  method  described  here,  there  are  no 
adaitional  parameters  other  than  those  determined  in  the  bond 
orb  tal  model. 
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FIGURE  CAPTION 


Fig.  1 The  normalized  temperature  rise  N(0,0,W)  at  the 
surface  (z  = 0)  and  beam  center  (r  = 0)  for  a laser 
oeam  of  shape  exp(-r2/w2'  which  attenuates  as 
exp(-az)  plotted  versus  W = aw.  The  normalization 
is  such  that  N(  oo  ) =1,  that  is,  N(W)  is  the 
ratio  in  temperature  rise  when  the  heat  is  depo- 
sited in  the  finite  attenuation  length  1/a  to  the 
rise  that  would  be  produced  if  the  same  heat  were 
absorbed  in  an  infinitesimal  layer  at  the  surface. 
The  upper  (solid)  curve  corresponds  to  the  horizon- 
tal scale  shown.  The  lower  (dashed)  curve  is  the 
upper  curve  replotted  with  the  horizontal  scale 
expanded  by  a factor  of  5 to  facilitate  reading 
points  in  the  small  w region. 
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